冷 热 油 交 蔡 输送 热力 快速 计算 研究 


韩 东 旭 , 王 月 ， 邵 倩 倩 ， 宇 波 

(中 国 石 油 大 学 机 械 与 储 运 工程 学 院 ， 城 市 油气 输 配 技术 北京 市 重点 实验 室 ， 北 京 102249) 
摘 要 : 本文 将 基于 贴 体 坐 标的 非 稳 态 导热 POD 低 阶 模型 应 用 于 冷 热 油 交替 输送 热力 快速 计算 中 ,并且 详细 介绍 了 了 POD 低 
阶 模型 的 实施 方案 。 在 低 阶 模型 的 求解 过 程 中 针对 管道 热力 影响 区 域内 各 部 分 导热 系数 相差 达到 300 倍 左右 ， 文 中 采用 更 符合 
基 函 数 空间 导数 物理 意义 的 离散 格式 代替 传统 的 中 心 差分 格式 来 对 其 进行 离散 。 在 本 文 的 算 例 中 POD 低 阶 模型 平均 计算 误差 
不 超过 0.35*C， 计 算 速 度 为 传统 的 有 限 容积 法 的 260 倍 左右 。 
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Abstract: A proper orthogonal decomposition (POD) Galerkin reduced order model based on body-fitted coordinate 
ls applied to the fast calculation of the thermal process during cold and hot oil alternating transportation, and the 
implementation process is given in detail. Considering the uneven distributed heat conductivity coefficient in the 
thermodynamically involved region of the pipeline (the ratio of heat conductivity coefficients is as large as 300 
around), a discretization Scheme better fitting the physical meaning of the derivative of the basic function Space is 
employed, instead of the traditional used central difference scheme. Trough all the cases devised, the averaged 
computation error by POD reduced order model is below 0.35°C, while the speed can reach 260 times of the 
traditional finite volume method. 
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calculation 
0 前 言 热力 规律 及 其 停 输 再 启动 问题 。 在 上 述 数值 研究 中 ， 
瑟 


存在 一 个 普遍 问题 即 计算 速度 较 慢 。 在 对 热 油管 ; 
冷 热 油 交替 输送 属国 际 前 沿 原油 输送 技术 ， 即 的 优化 设计 和 风险 评估 中 往往 需要 对 大 量 的 不 同 工 
对 两 种 或 多 种 原油 根据 其 物性 尤其 是 粘 温 关系 的 不 况 进行 计算 ， 这 就 会 耗费 大 量 的 时 间 ， 甚 至 使 得 计 
同 采 用 不 同 的 温度 进行 输送 ， 在 保证 安全 输送 的 前 算 不 能 完成 或 者 丧失 了 计算 结果 的 时 效 性 ， 所 以 必 
提 下 减少 加 热量 达到 节能 的 目的 。 但 是 这 就 会 使 整 须 寻 求 快速 的 计算 方法 。POD 低 阶 模 型 作为 一 种 快 
条 管线 始终 处 于 热力 不 稳定 状态 ， 为 了 保证 安全 输 速 计算 方法 被 广泛 应 用 于 各 工程 领域 中 1。 
送 需 要 对 其 热力 规律 进行 深入 研究 。 在 输油管 道 的 热力 规律 计算 方面 文献 [10] 首 次 
崔 秀 国门 在 其 博士 论文 中 从 理论 上 深入 研究 ] 将 POD 低 阶 模型 引入 到 了 热 油管 道 的 计算 中 ， 采 用 
这 一 输 油 方案 , 姚 峰 外 在 大 量 数值 模拟 的 基础 上 为 中 RO AAO dri 快速 计算 了 埋 地 热 
石油 西部 管道 公司 制定 了 整套 输送 工艺 , 宇 波 等 站 管道 的 冷 热 油 交 蔡 输 送 和 原油 投产 ; 文献 [11] 采 用 
EM 步 研究 了 冷 热 油 交 蔡 输 送 管道 的 。 基于 贴 体 坐 标的 稳 态 导热 POD 低 阶 模型 ， 通 过 贴 
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体 坐 标的 计算 平面 建立 起 不 同 的 埋 深 和 管 径 的 管 
横 截 面 的 关系 ， 通 过 一 次 取样 即 可 快速 测 不 同 
和 埋 深 的 管道 横 截面 稳 态 温度 场 。 

但 文献 [10] 和 [11] 均 未 考虑 埋 地 管道 的 结 蜡 层 、 
钢管 壁 和 防腐 层 ， 由 于 其 导热 系数 相差 巨大 ， 严 习 
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题 ; (4) 管内 结 蜡 层 厚度 在 管道 周 向 和 轴 向 均 不 发 
生变 化 ; (5) 管道 的 热力 影响 区 为 10m， 即 图 2 中 
L=10m，H=10m。 

结 蜡 层 、 钢 管 壁 、 防 腐 层 、 土 壤 的 传 热 控制 方 


a 


影响 POD 低 阶 模型 的 计算 精度 。 文 献 [10] 不 能 天 
同 埋 深 和 管 径 的 热 油管 道统 一 快速 计算 ， 文献 [11] 仅 
仅 考 虑 了 管道 横 截 面 的 稳 态 温度 场 。 

针对 以 上 问题 ， 本 研究 采用 基于 贴 体 坐 标 非 稳 
态 导热 POD 低 阶 模型 ， 采 用 特殊 格式 离散 基 函 数 的 
则 导数 ， 成 功 快 速 计算 了 不 同 埋 深 和 管 径 的 热 油 
管道 在 冷 热 油 交 蔡 输 送 时 的 热力 规律 。 


1 冷 热 油 交 丛 输 送 物理 模型 


图 1 所 示 为 埋 地 冷 热 油 交 蔡 输 送 管道 示意 图 ， 
冷 油 和 热 油 交 蔡 从 热 站 流入 管道 ， 油 流 与 结 蜡 层 内 
侧 对 流 换 热 ， 结 蜡 层 、 管 壁 、 防 腐 层 以 及 土壤 以 热 
传导 的 形式 换 热 ， 并 经 地 面 与 大 气 换 热 。 图 2 为 埋 
地 热 油管 道 的 热力 影响 区 。 
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图 1 埋 地 冷 热 油 交 蔡 输 送 管道 示意 
Fig. 1 Sketch of buried pipeline for cold and hot oil alternating 
transportation 


i 
管道 x 方向 热 
力 影响 区 2L 
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图 2 管道 横 截 面 示意 图 
Fig. 2 Schematic diagram of computational domain 


为 了 方便 计算 ， 在 冷 热 油 交 蔡 输送 油 温 沿 程 分 
布 的 计算 过 程 中 ， 做 出 以 下 假设 :(1) 管内 原油 温 
度 在 某 一 截面 上 是 均匀 的 ;(2) 管道 周围 的 土壤 为 
各 向 同性 的 均匀 介质 ; (3) 忽略 土壤 区 域 轴 向 温 降 ， 
将 三 维 非 稳 态 导热 问题 简化 为 二 维 非 稳 态 导热 问 


程 : 
oT 0 oT、 6 oT 
C 三 4 4 K =12,3,4 
Pr k dz Be k Bx By k Ea (1) 
油 流 控制 方程 
连续 性 方程 : 
0 0 
一 (OA)+ 一 (OVA) =0 
pe A PN CO) 
油 流 的 换 热 方程 ; 
d9 9 dp fVi 49 
C 
”dr 2 2D pD G3) 


油 流 与 结 蜡 层 之 间 的 匹配 条 件 : 
oT 
We 三 下 公 - 生 ) (4) 
时 中 了 ,表示 结 蜡 层 和 油 流 接触 面 的 温度 ，9 表示 热 
油 温度 , h, 是 管内 油 温 对 结 蜡 层 的 强制 对 流 换 热 , h 
为 管内 油 温 的 函数 。 
管道 热 影响 区 域 边界 条 件 : 


x=0,-(H,-R+6)<y<0 时 ， 


芝 


党 和 三 0 当 
Ox 


x=0,-H<y<-(H,+R-6) 时 , 和 二 -0 当 
X 


oT si 
y=0 时 ? hh Toroma) 3 当 X| = 工 时 》 
0 ge 
Ox 


其 中 ，h, 为 大 气 对 土壤 的 换 热 系数 ， 其 值 由 实 
际 测 量 确定 ， 工 为 土壤 恒温 层 温度 其 值 由 实际 情况 
确定 ， 也 为 大 气温 度 ， 在 本 文中 采用 式 〈5) 来 计算 


， ,元 
I 二 In ti)+t+at, +b (5) 


式 (5) 中 ，8 单 位 为 小 时 ，Ssin(176) 为 一 天 内 
气温 变化 ; 5 单位 为 天 ，ab +b 为 气温 随 日 期 的 变化 。 
2 直接 数值 模拟 


在 冷 热 油 交替 输送 的 热力 计算 中 ， 主 要 需要 求 
解 方程 (1) 和 (3)， 下 面 对 这 两 个 方程 的 直接 数值 
求解 做 简要 介绍 。 

求解 方程 1) 主要 目的 在 于 求解 任 一 管道 横 截 
面 处 管道 热力 影响 区 域内 的 温度 场 ， 根 据 对 称 性 为 
了 减少 计算 量 只 计算 图 2 所 示 区 域 的 一 半 。 本 文中 


采用 基于 贴 体 坐 标的 FVM 计算 热力 影响 区 域内 的 温 
度 场 ， 当 求 解 区 域 转换 到 了 计算 平面 后 ， 控 制 方程 
相应 的 转换 到 了 计算 平面 中 ， 控制 方程 公式 (1) 转 
换 成 公式 (6) 
OPCT) _ 0 74 


加 5 人 (oT — PT)]1 


0 (6) 
07 
在 FVM 计算 过 程 中 , 采用 的 贴 体 网 格 的 总 网 格 
数目 为 3761， 不 同 物性 界面 的 导热 系数 采用 调和 平 
均 进 行 插值 ， 采 用 Gauss-Seidel 迭代 ， 连 代 收 敛 标准 
为 相 邻 两 个 迭代 层 的 对 应 控制 点 最 大 误差 小 于 10”。 

利用 原油 与 周围 环境 热量 平衡 关系 实现 管内 油 
流传 热 与 管 外 土壤 传 热 的 耦合 ， 采 用 特征 线 法 求解 
方程 〈2) 并 获得 不 同时 间 沿 程 油 温 分 布 。 


3 采用 POD 低 阶 模型 快速 计算 


冷 热 油 交 蔡 输送 的 热力 规律 计算 ， 其 计算 量 主 
要 集中 在 管道 横 截 面 上 温度 场 的 计算 即 方程 (6) 的 
求解 ， 该 部 分 的 计算 消耗 时 间 占 据 整个 计算 时 间 的 
99.5% 以 上 。 为 此 本 文 提 出 采用 基于 贴 体 网 格 的 POD 
低 阶 模型 来 快速 求解 该 方程 ， 从 而 实现 对 不 同 埋 深 
与 管 径 的 冷 热 油 交 蔡 输送 热力 规律 的 快速 计算 。 
由 POD 定义 可 知 ， 稳 态 时 计算 平面 上 的 温度 场 
可 写成 一 组 基 函 数 和 对 应 的 谱系 数 的 线性 县 加 的 形 
式 ， 其 中 M 表示 所 有 基 函 数 的 个 数 。 


T(6,7) = 2 arg ($7) (7) 


1 
0 PT)+JS k=1,2,3,4 


基于 方程 (6), 残 差 R 可 以 定义 为 ， 
O(pc,T) 6 .4 


A) 
OO 4 
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将 方程 (89) 向 基 函 数 的 空间 投影 ， 经 过 一 系列 推导 
可 得 〈 详 见 文献 [12]) ， 


M da， M M 
Ve i ee 
2 时 =1 =1 


中 VcaqsWdE+ Vrq'gbdnt+r, 


(9) 


E 0 从 0 从 i 0 从 0 从 
RE 


F=-(S,0,), ©;,=%:J 

将 初始 条 件 和 边界 条 件 引 入 到 方程 (9) 中 可 求解 
得 到 谱系 数 ， 详 见 文 献 [12]。 
由 于 管道 热力 影响 区 域 中 ， 结 蜡 层 、 钢 管 壁 、 
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防腐 层 和 土壤 层 的 导热 系数 (本 文中 分 别 取 0.15 
w/(m'C)、 50 (wmC)、0.15 (w/m'C) 和 0.85 (w/m'C)) 
相差 巨大 ， 尤 其 是 钢管 壁 同 结 蜡 层 和 防腐 层 的 导热 
系数 相差 300 倍 左右 。 在 这 种 情况 下 ， 采 用 传统 的 
中 心 差 分 格式 来 对 其 进行 离散 已 经 不 能 满足 精度 要 


求 ， 本 文采 用 与 基 函 数 空间 导数 物理 意义 相 适 应 的 
如 下 离散 格式 ， 
(从 ) 办 +Br /Bo 从 Gp + By /Bp 办 
$7? 1+B,/B, 1+ By / Bp (10) 
(办 ) 二 Gp + By / Bp,py Gp + Bs / Bp,gs 
”? 1+B,/B,, 1+ B, / B,, (11) 
A Aa Aa 
By = iBs = i B = 9» 
其 中 JV7l, JVY | J Va lw 
Aa A 4Cw 
Brz = :Ba = ,Bp, = 
” JVc | J Va | JVY|, 


该 离散 格式 根据 调和 平均 推导 得 到 ， 可 更 好 地 
保证 能 量 守 恒 ， 因 此 采用 该 格式 离散 基 函 数 空间 导 
数 的 低 阶 模型 具有 较 高 的 精度 《〈 详 见 文献 [12] )。 


4 算 例 分 析 


在 冷 热 油 交 蔡 输送 管道 的 设计 ， 以 及 运营 方案 
的 优化 中 需要 对 大 量 的 工 况 进行 热力 数值 模拟 ， 直 
接 数 值 模拟 通常 耗费 时 间 较 多 ， 下 面 就 以 一 个 实际 
问题 为 例 介 绍 采 用 POD 低 阶 模型 是 如 何 快 速 求解 这 
一 类 问题 的 。 

两 条 管线 pipe 1 和 pipe 2， 为 两 条 已 经 达到 热 
力 平衡 的 稳 态 运行 输油管 线 ， 两 条 管线 埋 深 、 管 径 
甚至 钢管 壁 厚 都 有 很 大 的 差异 ， 各 项 运行 参数 见 表 
1。 根 据 工程 实际 的 要 求 ， 需 将 其 中 一 条 改 为 冷 热 油 
交 蔡 输送 工艺 ， 由 于 输送 方案 较 多 ， 计 算 量 较 大 所 
以 采用 基于 贴 体 坐标 低 阶 模型 来 快速 计算 不 同方 案 
的 热力 规律 ， 下 文 详 细 介 绍 了 其 实施 步骤 。 
表 1 管线 基本 参数 


Table 1 Pipeline parameters 


4 


管线 站 间距 埋 深 管 径 辟 睹 进 站 温度 
/km /m /m /m /'C 
pipe 1 108 13 0.711 0.010 50 
pipe 2 108 0.9 0.457 0.006 30 


为 了 说 明 低 阶 模 型 的 精度 和 计算 速度 ， 选 取 众 


多 方案 中 具有 代表 性 的 两 个 来 对 低 阶 模型 进行 测 
试 : 方案 1 更 改 pipe 1 为 30C 和 60C 交 蔡 输 送 ，10 
小 时 交 蔡 一 次 ， 流 速 为 1m/s， 在 夏季 进行 更 改 ， 即 
式 (5) 中 b=20,a= 0.4; 方案 2 更 改 pipe 2 为 20C 和 
50C 交 蔡 输 送 ，20 小 时 交 蔡 一 次 ， 流 速 为 2ms， 在 
冬季 更 改 ， 即 式 (5) 中 b=-20,a= -0.4。 


首先 进行 取样 。 低 阶 模型 的 主要 目的 是 计算 冷 
热 油 交 蔡 输 送 时 管道 横 截 面 的 温度 场 ， 所 以 取样 目 
的 是 为 了 获得 管道 横 截面 的 温度 场 信息 ， 各 参数 的 格 采 用 第 


取样 值 见 表 2 和 表 3, 取样 时 各 样本 从 初 场 一 直 计 算 
到 稳 态 ， 时 间 步 长 为 1200s， 在 初始 阶段 温度 场 变 化 
剧烈 时 加 大 取样 密度 而 在 温度 场 变 化 缓慢 的 后 期 则 
应 减 小 取样 密度 。 由 于 取样 只 是 针对 一 个 管道 的 横 
截面 进行 数值 计算 而 不 是 对 一 条 管道 的 所 有 截面 进 
行 计算 ， 所 以 取样 付出 的 时 间 代 价 较 小 。 
其 次 求 取 基 函数 。 在 获得 样本 和 天 阵 之 后 采用 快 
照 法 计算 基 函 数 ,前 8 组 基 函 数 包 含 了 基 函 数 总 能 量 
的 99.96%。 由 于 几何 参数 变化 的 影响 ， 导 致 基 函 数 


行 时 ， 


而 此 时 管内 进口 原 ; 
只 需 15 个 小 时 就 可 以 将 整个 管道 内 的 
毕 。 无 论 是 对 于 方案 1 还 是 方案 2, POD 低 阶 模型 都 
表现 出 了 较 高 的 精度 
通常 情况 下 ， 我 们 最 关心 的 是 


沿 程 


进行 计算 ， 二 者 采用 的 时 间 步 长 都 是 1200s， 计 入 
2 节 所 述 网 格 ，FVM 的 迭代 方法 和 收敛 标 
准 也 如 前 所 述 。 
图 3 给 出 了 方案 1 和 方案 2 的 新 工艺 在 初期 运 
具有 代表 性 的 时 刻 原 


中 。 采 用 FYM 和 POD 低 阶 模型 分 别 对 这 两 种 方案 


网 


温度 分 布 图 。 方 


案 1 需要 30 个 小 时 才能 将 原来 管 中 的 油 置换 完全 ， 


已 经 发 生 了 两 次 交 蔡 ;方案 2 


y 
个 


置换 完 


管道 内 原油 的 进 


中 可 以 看 出 


中 一 些 较 重要 的 局 部 信息 存在 于 后 面 能 量 贡 献 较 低 。 站 温度 ， 为 此 图 4 分 别 给 出 了 pipe 1 和 pipe 2 在 30 
的 基 函 数 中 ， 所 以 在 POD 低 阶 模型 计算 中 , 为 了 准 ” 天 内 进 站 温度 随时 间 变 化 的 曲线 。 从 恬 
确 计 算 管 道 横 界面 的 温度 场 ， 采 用 了 前 26 个 基 函 数 ” 方案 1 和 方案 2 中 进 站 油 温 周期 性 变化 ， 周 期 分 别 
来 对 其 进行 预测 。 为 10h 和 20h， 与 出 站 油 温 的 周期 性 变化 一 致 
最 后 ， 将 低 阶 模型 应 用 于 前 述 的 两 个 测试 方案 
表 2 冷 热 油 交替 输送 管线 横 截 面 取 样 
Table 2 Samples on the cross section of the pipeline 
编 油 温 油管 换 热 系 数 空气 温度 大 气 换 热 系数 初 几何 
号 /°C /W/(m”.C) /°C /W/(m”.C) 场 形状 
1 20 25 -5 10 空气 温度 为 0C， 人 恒温 
2 40 55 0 25 层 温 度 为 SC 的 稳 态 温 。 ”Geol1, 见 表 3 
3 60 90 15 25 度 场 
4 20 25 -15 10 空气 温度 为 0C， 人 恒温 
5 40 55 0 10 层 温 度 为 -5C 的 稳 态 温 。 Geo2, 见 表 3 
6 60 90 5 25 度 场 
4 20 25 15 25 空气 温度 为 10C， 恒温 
8 40 55 0 层 温 度 为 -5'C 的 稳 态 温 。 Geo3, 见 表 3 
9 60 90 -15 25 度 场 
表 3 冷 热 油 交替 输送 取样 表 
Table 3 Samples for alternating transportation of cold and hot oil 
埋 深 (m) 管 径 (m) 壁 厚 (m) 结 蜡 层 (m) 防腐 层 (m) 
Geol 1.5 0.8 0.007 
Geo2 1 0.5 0.012 0.022 0.004 
Geo3 13 0.6 0.09 


Oil emperature (C) 
Oil tempreure ('C) 


mpreture (C) 


(a) 方案 1 
图 3 投产 油 温 沿 程 分 布 图 
Fig. 3 Oil temperature distribution along pipeline 
during the commissioning process 


(b) 方案 2 


(a) 方案 1 


(b) 方案 2 


图 4 ; 


站 温度 随时 间 变 化 图 


Fig. 4 Temperature at the inlet of the station 


从 图 4 中 可 以 看 出 POD 计算 结果 与 FVM 计算 
结果 吻合 良好 ， 两 种 方案 POD 计算 最 大 误差 不 超过 
0.9C， 平 均 误 差分 别 为 0.35'C 和 0.30'C， 这 个 误差 
水 平 在 输 油 行业 是 完全 可 以 接受 的 。 

为 了 更 加 形象 的 说 明 POD 低 阶 模型 的 计算 精 
度 ， 图 5 (a) 和 图 5(b) 分 别 给 出 了 pipel 末端 进 站 处 
在 冷 热 油 交 蔡 输送 的 第 5 天 和 pipe2 末端 进 站 处 第 


(a 方案 1 第 5 天 


(b) 方 案 2 第 30 
图 5 进 站 处 管道 横 截 面 温度 分 布 图 


Fig. $ Temperature distribution on the cross section at the inlet 


为 了 说 明 相 比 于 直接 数值 模拟 的 FVM 方法 ， 
POD 低 阶 模型 具有 较 高 的 速度 优势 ， 本 文 测算 了 在 
处 理 器 为 Intel i7 CPU @3.50GHz 的 微型 计算 机 上 ， 
两 种 方法 计算 所 耗费 的 时 间 : 对 于 方案 1，FVM 耗 
时 3900s， 采 用 POD 耗 时 14.8s，POD 计算 速度 为 
FVM 263 倍 ; 对 于 方案 2，FVM 耗 时 2220s， 采 用 
POD 耗 时 8.7s,，POD 计算 速度 为 FVM 255 倍 。 需要 
注意 的 是 , 方案 2 的 计算 时 间 少 于 方案 1， 是 由 于 方 
案 2 流速 快 ， 需 要 计算 温度 场 的 截面 较 少 。 假 如 工 
程 中 有 200 种 方案 (各 参数 不 同 取 值 的 排列 组 合 ， 
通常 这 个 数量 较 大 ) 可 供 选 择 ， 如 果 采 用 FVM 则 大 
约 需要 耗费 7 天 的 计算 时 间 ， 这 在 工程 上 是 几乎 不 
能 承受 的 ， 而 POD 低 阶 模型 则 只 需要 40 分 钟 ， 这 
为 工程 上 优化 设计 和 运营 都 带 来 了 极 大 的 方便 


5 结论 


本 文采 用 基于 贴 体 坐标 的 非 稳 态 导 热 POD 低 阶 
模型 ， 基 函数 空间 导数 采用 更 符合 其 物理 意义 的 格 
式 离 散 ， 实 现 了 对 不 同 管 径 和 埋 深 管道 冷 热 油 交 蔡 
输送 热力 规律 的 快速 计算 。 测 试 算 例 表明 低 阶 模型 
有 较 高 的 精度 和 健壮 性 ， 完 全 可 以 满足 工程 需求 ， 
且 其 计算 速度 在 本 文 算 例 中 为 FVM 方法 的 260 倍 左 
右 。 该 快速 计算 方法 的 应 用 会 大 大 加 快 工程 中 冷 热 
油 交 蔡 输 送 管道 的 优化 设计 、 运 营 以 及 风险 评估 的 
速度 。 
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